Landau-Zener tunneling with many-body quantum effects in crystals of molecular 

magnets 
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We present a quantum interpretation of the heights in hysteresis of Fes molecule at lower tem- 
peratures by treating the crystal as an Ising spin system with the dipolar interaction between spins. 
Then we apply it to two limit cases : rapid and adiabatic regions. Our theoretical analysis is in 
agreement with the experimental observation in these regions, which indicates that the steps in 
hysteresis loops of magnetization of Fes at lower temperatures show a pure quantum process. 
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INTRODUCTION 

Crystals of molecular magnets, such as Fes and Mni2, 
have attracted much attention for their connection to 
many macroscopic quantum phenomena QQiSH- They 
may also have important applications in magnetic mem- 
ory and quantum computing |4j,|5j. The earliest and most 
spectacular observation on such a system is the quan- 
tum steps in the hysteresis loops of magnetization at low 
temperatures 0, y[ . These quantum steps are the man- 
ifestations of macroscopic quantum tunneling, resulting 
from the tunneling between different spin states of large 
molecular spins (S = 10 for both Fes an d M1112, S = 9/2 
for M114). Fes is particular interesting because the steps 
in the hysteresis will become temperature independent 
below 0.36-ftT, which shows a pure tunneling process 0. 
This tunneling phenomenon is complicated by the in- 
teraction between spins and other environmental effects. 
Intensive efforts have been devoted to explain the step 
features via different approaches 0, 0, U The modi- 
fications of other environmental effects have also been 
studied, such as the nuclear spin effects [9j, [l(J, [llj . 

More recently, Liu, et al presented a successful the- 
ory on the height of quantum steps in the hysteresis loop 
when the temperature is low enough that the thermal 
effects can be neglected 0. By treating Fe 8 crystal as 
a system of Ising spins sitting at each site of the lattice, 
the step heights measured in experiment were success- 
ful reproduced by directly solving an evolution equation 
described by the dipolar fields distribution. The results 
have been compared to Landau-Zener (LZ) model |l3| . 
which has been used to extract the tunnel splitting A 
of a single molecular spin from step heights 0, ^| . In 
their simulation, the dipolar interaction between spins is 
treated by a mean-field theory and the flipping of each 
spin is independently. 

In this paper, we also model the Fes crystal as a system 



of Ising spins sitting at each site of the lattice and taking 
into account the dipolar interaction between spins. But 
different form Ref. we treat the system as a quan- 
tum many-body system and formally give a formula to 
evaluate magnetization. Although it can not be used to 
calculate magnetization in most cases for the algorithm 
reason (the time of calculating is increasing exponentially 
with the number of sites), we apply it to two limit cases: 
rapid and adiabatic cases. Through the pure quantum 
approach, our theory successfully interpret the quantum 
step heights in these two regions. As an application of our 
theory, we show that the tunnel splitting A e measured 
with the LZ model 0, ^| is proportional to the true 
tunnel splitting A of a single molecular with a geometry 
factor (depending on shape and lattice structure). This 
result has also been obtained in Ref. [12j. In the adia- 
batic limit, we show the measured tunnel splitting A e is 
dependent with the sweeping rates with the power law: 
A e ~ a 1 / 2 . This prediction agrees with the experimental 
observations. 



MODEL 

We model the Fe 8 crystal as spin lattices with the real- 
istic constants 10.52 : 14.05 : 15.00 (A) and angles 89.9° : 
109.6° : 109.3° between the axes, which is a triclinic lat- 
tice; the shortest axis a as the easy axis (actually there 
is an angle of about 8° between them, but it does not af- 
fect significantly the results). As in Ref. we focus on 
one step for simplicity, that is, the tunneling between the 
two lowest levels (S z = ±10). The effective Hamiltonian 
operating in the subspace is j^l 

W = - E MBStoHa? - \ E Vij of a® + ± E A(J * l) ■ 

i ij i 

(1) 
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The first term describes the Zeeman energy, and H is 
the external field applied in the direction of the easy 
axis. The second term is the spins interaction with 
dipolar potential Vij — E d (3cos 2 9 — l)Q /rfj, E d = 
^■(5/is5) 2 /f2o, where is the displacement between 
the spins, 9 is the angle between f\j and the easy axis, 
and f^o is the unit cell volume. The last term describes 
tunneling and A is the tunnel splitting. a z and a x are 
Pauli matrices, and {i}, {j} label molecular sites. 

Because A/E d ~ 1CT 6 (see Ref. [l5|), the last term of 
Eq. can be regarded as a small perturbation, denoted 

as W = \ J2i ■ We can know the eigenstates of 
the unperturbed system are (p^ = | s±, s% , • • • , Spj) (a* — 
1, • • • , 2^), in which — ±1 corresponding to the spin 
on site i up and down respectively. The eigenvalue for 

<t>„. is 



(X (J, 

13 S i S j ■ 



(2) 



where P_(jy-i) is the sum of probability for all the en- 
ergy levels with i spins up. For example, assuming there 
are m levels of one-spin up states, which are denoted as 



E' 
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(1 = 1,2,- 



The energy structure is shown 



in Fig.l where only the initial level E_n and one-spin up 
levels are plotted. The energy gap between E l _, N _y and 
E-n is 5[A. Then we can obtain the probability of the 
initial level, 



P-N = 



l[Ps[ 



= P. 2a 



(4) 



where c = Y^iLi(^l) 2 - I* 1 analog, we can evaluate the 
probability for any level involved in principle. It can be 
formally expressed as 



(N-l) 



E 



Y[0--ps,)]Ip5* 

P=l i 



(5) 



It is easy to see that the energy levels would be degener- 
ate when some sites are geometric equivalent. 

On the other hand, because operating on M will 
make the spin on the z-th site flip, the perturbation term, 
i.e., the off-diagonal element, 7i M „ = (frf/H^v = 4>^W(j) v , 



is not zero if and only if s. 



{i = 1.2 j 1:./ ' 



l,j + 2, ...,N) and Sj ^ Sj, where j could be any site. 
At this time = A/2. 

From the perturbation theory |lf| , we can know, if off- 
diagonal element of two states is nonzero, the correspond- 
ing energy levels must have an avoided crossing with a 
gap proportional to the off-diagonal element A/2. The 
gap of the avoided crossing is determined by the degen- 
erate properties of the levels. For example, if both levels 
are non-degenerate, the gap is A; but if one of them is 
two- fold degenerate, the gap is v2A, and so on. Higher 
order perturbations are much small, so they can be re- 
garded as crossings. For example, for the second order, 
the gap is about A 2 /E d ~ 10 _6 A, so it can be treated 
as a crossing. 



in which we have ignored some subscripts for conve- 
nience. 
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FIG. 1: The energy structure of the initial level E-m and one- 
spin up levels. The energy gap between E l _^ N _ 1 - ) and E-n is 

S[A. 



MAGNETIZATION AND TUNNEL SPLITTING 
Result of our model 

Supposing the crystal has N Fes molecules, it is ini- 
tially on the state | — 1, — 1, • ■ • , — 1) in a large negative 
field, then sweep the field with a constant rate a to the 
positive. Over an avoided crossing, the spin involved will 

7r((SA) 2 

flip with the probability 1 — p$ where p$ = e 2° if the 
gap is S A. The magnetization can be formally expressed 
as 

M = -NSP- N - (JV-2)5P_ (JV _ 1) + 

+ ..._(JV-2*)SP_ (JV _ <) + ..., (3) 



There should be a term in Eq. corresponding to the 
adiabatic path: starting from the initial level and keep- 
ing on the continuous branch at every avoided crossing 
encountered. The probability of the adiabatic path can 
be expressed as 

P ad = (l-p sl )(l-p sl )---(l-p 5 i K ), (6) 

where (1 — p^i) is the probability of spin flipping over 
each avoided crossing and K is just the total number of 
the flipped spins along this path. 

Indeed, can not be calculated for nowaday com- 
puter when the number of spins is large. But we can 
evaluate it for limit cases. 
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For the high sweeping rate limit, P = e - " 2 ^ — » 1, so 
(1 — P) is an infinitesimal. Form Eq. (J5J, one knows 
that P_(jy-i) is as the same order as (1 — P) 1 . Hence, 
to the first order of (1 — P), the magnetization can be 
approximately expressed as 

M « -NSP-n - (N - 2)5P_ ( jv_i). (7) 

Of course, the total probability must be conserved, i.e., 



P_ 



N 



P_ 



(JV-1) 



P_ 



(N-i) 



= 1. 



(8) 



So, to the first order of (1 — P), we have P_n + P_ 



(JV-1) 



1, i.e., P_ 



(JV-l) 



1 - P_ 



Substituting (0J into the 



above formula, we get 



M « -iVS* + 2&S'- 



2a 



(9) 



In Refs. U,llJ,ll5|i Wernsdorfer and coworkers extract 
the tunnel splitting A of a single molecular spin from the 
magnetization by employing LZ model 01 ■ Based on LZ 
model, the measured tunnel splitting A e was calculated 
by the following formula [IsT ]. 



(10) 



in which M s = NS. 

Substituting the theoretical predication of magnetiza- 
tion Q into the above formula, we obtain 




(11) 



in which C = 

This result shows that for the rapid sweeping rate limit, 
the measured tunnel splitting A e is a constant, which 
consists with the experiment observation (lif. Eq. i|ll|) 
also implies that A e is not the true tunnel splitting A, 
but proportion to A with a factor C which is dependent 
only on the geometry of the sample: its shape and lattice 
structure. This consists with the result presented in Ref. 

If the sweeping rate is very small (adiabatic limit), 



p S j -> 0, so P 

Pad — ► 1, SO M 

we obtain 



-(N-i) ~~ > except for the adiabatic term 
M ad = -{N - 2K)S, then from JUJ 



ka 1/2 , (12) 




Comparing our results with experiments 

For the high sweeping rate limit, as shown in Eq. I jllj l. 
the measured tunnel splitting is a constant, which con- 
sists with experimental observation. For the adiabatic 
region, we can calculate the measured tunnel splitting in 
the adiabatic limit since the number of levels involved 
in adiabatic path is proportion to the number of spins 
N. We can find the adiabatic path by following the adi- 
abatic process: starting from the initial level and keep- 
ing the state on the continuous branch at every avoided 
crossing encountered. In Fig. 2, we plot the adiabatic 
magnetization M a d of n x n x n lattice for different total 
number of spins. One can find that as the total number of 
spins is large enough, M a d becomes independent on the 
total number, and tends to a constant M a d/M s ~ —0.29, 
i.e., k = 0.528 . We also calculate M a d for the case of 
(axbxc) : 16 x 8 x 8 lattice (a is easy axis direction), and 
obtain M ad /M s = -0.37, i.e., k = 0.49. In Fig. 3, we 
compare the theoretical evaluation with the experiment 
observation 0] . It is shown that for three different Pe§ 
isotopes, the three curves of A e tend to merge together 
with the same tendency A e oc a 1 / 2 . This tendency con- 
sists with our theoretical prediction. We argue that the 
adiabatic evolution of the system is only determined by 
the levels structure, so we cannot read the information 
of the tunnel splitting from the adiabatic process. This 
feature can be found in Fig. 5 where three curves of 
different isotopes have the same tendency. 




FIG. 2: The adiabatic magnetization ofnxnxn Pes triclinic 
crystals for different number of spins. The solid line is guide 
for eyes. iV = nxnxn is the total number of spins. 



in which k 



■In 



1-M ad /M s 
2 



This shows that for 



the slow sweeping rates, the measured tunnel splitting 
is strongly dependent on the sweeping rate. In the adia- 
batic limit it shows a 1/2 power law of the function of the 
sweeping rate: A e = fca 1 / 2 . This feature is first revealed 
in this paper. 



CONCLUSION 

In summary, we have given a pure quantum interpreta- 
tion of the step heights in hysteresis loops of Pes molecule 
by treating the crystal as a system of Ising spins sitting 
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FIG. 3: A e of Fes crystals for small sweeping rates. The solid 
line is theory prediction with n x n x n, the dashed line is for 
2n x n x n, the others are the experimental data for different 
Fes isotopes [15]. 



at each site of the lattice with the dipolar interaction 
between spins. Our theoretical analysis is in agreement 
with the experimental observation in the rapid and adia- 
batic limits. For the rapid sweeping rates, we show that 
the measured tunnel splitting A e is a constant which is 
proportional to the tunnel splitting of the single molec- 
ular spin, i.e., A e = CA. The factor C depends on the 
sample geometry. But for the adiabatic limit, the magne- 
tization becomes independent on the sweeping rate, and 
tends to a constant. This feature leads to that the mea- 
sured tunnel splitting A e is strongly dependent on the 
sweeping rate a, and be a 1/2 power law of sweeping 
rate: A e ~ a 1 / 2 . 
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